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\Q ! ABSTRACT 

O 
O , 

CN . Core collapse of massive stars and binary neutron stars or black hole-neutron star binary 

mergers are likely progenitors of long and short duration gamma-ray bursts respectively. Neu- 
tronized material in the former and neutron star material in the latter are ejected by the central 
Q\ \ engine implying a neutron-rich jet outflow. A free neutron, however, beta decays to a proton, an 

_ -, ' electron (beta) and an anti-neutrino in about fifteen minutes in its rest frame. Sudden creation of 

CN ■ 

J> ■ a relativistic electron is accompanied by radiation with unique temporal and spectral signature. 
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We calculate here this radiation signature collectively emitted by all beta decay electrons from 



£T) . neutron-rich outflow. Detection of this signature may thus provide strong evidence for not only 



neutron but also for proton content in the relativistic gamma-ray burst jets. 

Subject headings: gamma rays: bursts — gamma rays: theory — ISM: jets and outflows — radiation 
mechanisms: non thermal 
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' 1. Introduction 

C3 ■ 

The constituents of relativistic gamma-ray burst (GRB) jets are not known completely as yet. The 
observed MeV energy 7-rays are most likely radiated by energetic electrons which are shock accelerated 
by internal collisions of successively ejected materials by the GRB central engine. In this widely accepted 
fireball shock model, protons along with electrons are needed to be present in the GRB jet to explain the 
observed rapid variability of 7-rays and the GRB afterglow. Free neutrons are speculated to be present in 
the GRB jet (Derishev, Kocharovsky & Kocharovsky 1999a) which may be a common feature (Beloborodov 
2003b) following the ejection of neutronized core material in case of a long burst which are most likely 
created by core collapse of massive stars, a mechanism similar to that of the supernovae of type Ib/c or type 
II. Recent observations of several short burst afterglows suggest a common relativistic jet feature in both the 
long and short bursts. In the latter case, neutron star material may also feed the jet with free neutrons. 

Observational signatures (both electromagnetic and neutrino) of a neutron component in the GRB 
jets have been discussed by many authors, in the early fireball evolution phase (Derishev, Kocharovsky 
& Kocharovsky 1999a; Bahcall & Meszaros 2000; Meszaros & Rees 2000; Razzaque & Meszaros 2006) 
and later in the internal shocks or afterglow phase (Derishev, Kocharovsky & Kocharovsky 1999b; Pruet 
& Dalai 2002; Beloborodov 2003a; Fan & Wei 2004; Fan, Zhang & Wei 2005; Dermer & Atoyan 2006). 
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However, most of these predictions depend heavily on the jet model parameters and interpretation of data 
may be difficult. 

In this paper we discuss a generic signature of neutron-rich jets, namely the radiation emitted by rel- 
ativistic electrons created by neutron beta decay: n — > pe~v e in the GRB jets. This radiation signature 
should be ubiquitiously present in the electromagnetic signal detectable from all GRBs if the jet contains 
free neutrons. We also find that the signal strength do not vary rapidly for reasonable ratios of the neutrons 
to protons in the jet and/or other jet parameters. Thus detection or non-detection of this signature may serve 
as a powerful discriminator between GRB jets with and without a neutron component. 

We discuss the neutron-rich jet models of long and short bursts in Section 2, beta decay electron energy 
and number distribution in Section 3 and radiation signature in Section 4. We discuss our results and 
detection prospects of beta decay radiation in Section 5. 



2. Neutron-rich jet models 

The typical parameters we use for long and short GRB jets are listed in Table 1 . The baryon loading 
parameter (also known as the dimensionless entropy) is defined as the total energy to mass flow ratio in the 
GRB jet: rj = L/Mc 2 . The total mass outflow rate in the jet, ignoring the electron mass (m e <C m p ) and 
negligible thermal energy, is M ~ Airr 2 c(l + ^,o)n' p m p . Here £ G = n' n /n' is the initial neutron to proton 
number density ratio in the jet comoving frame. The outflow starts at a radius R Q which is a few times the 
Schwarzschild radius r g = IGM^/c 2 of a solar mass (few solar mass) black hole for a short (long) GRB. 

Electrons, protons and neutrons are coupled to thermal radiation in the expanding jet outflow as long 
as the Compton scattering time scales t' Th ~ (npCTThc) -1 and elastic n-p scattering time scale t' np ~ 
{n' p (Jnpc)~ l are shorter than the plasma expansion time t' cxp ~ r/cT(r) in the jet comoving frame (see, 
e.g., Rossi, Beloborodov & Rees (2005) for a detailed description of the neutron-rich jet dynamics). Here, 
T(r) oc r/R Q is the bulk Lorentz factor of the expanding jet outflow and the elastic n-p scattering cross- 
section a np dxh/20, the Thomson cross-section. The final value of F, after the expansion phase is over, 
depends on the value of rj. For sufficiently high value of 77, the neutrons may decouple from the outflow 
(Derishev, Kocharovsky & Kocharovsky 1999a; Bahcall & Meszaros 2000). The corresponding critical n 
value for n-p decoupling, from the condition t' np = t' cxp , is 

_ f La np 1 1/4 ^ |264 (L 52 / J R , 7 ) 1 /4(i + ^-1/4 

Vnp ^[4nR m p cHl+Z )\ ~ \l48 (L6o/^ 6 ) 1/4 (l +^i)" 1/4 
Here we have used £ G = 10£ Oj i. The final bulk Lorentz factor of the neutrons in the jet outflow is 

„ I Rnp/ Ro j V > Vnp 

1 nj = < (2) 

[T p j far); n < r) np 

where R np is the n-p decoupling radius where the nuclear scattering optical depth r' np ~ n' p a np R np /T{R np ) - 
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1 and is given by 

R (VnA 1/3 f 2.5 x 10 9 (L 52 Ry m . 6 )V3(i + e,i)- 1/3 cm 

V UJ ~\l.2xl0 8 (L 50 i2V^. 5 ) 1/3 (l+^i)- 1/3 cm 

Note that for 77 > r/ np the final Lorentz factor for the electron, proton and radiation outflow F p j > F n j 
[see, e.g., Rossi, Beloborodov & Rees (2005)]. In any case F n j < i], always. 

The kinetic luminosity of the freely coasting neutron outflow in the GRB jet is given by 

V 1 + £o 

Initially, at time t ~ 0, the neutron volume number density in the outflow, in the jet comoving frame, is 

Here (3 n ~ 1. However, this decreases exponentially with time as neutrons decay (beta decay) which we 
discuss next. 



3. Beta decay 

A free neutron has a mean lifetime Tp = 886.7 s in its rest frame. This corresponds to a mean decay 
radius Rp = crpT n j for a neutron at rest in the Lorentz boosted outflow. The volume number density 
of beta decay electrons (or protons or neutrinos) in the outflow evolves with time as n' e (t) = n' n o (l — 
exp[—r(t)/Rp]) and the total number of beta decay electrons at a time t is N' e (t) = Airr 2 (i)n' e (t)ct with 
no energy injection at later time t > 0. We denote the variables in the frame comoving with the neutron 
outflow with primes from now on. The relationship between the observed time and neutron outflow radius 
isr(t)=2ctTl J /(l + z). 

For our purpose we are interested in the total number of freshly created electrons at a time t which we 
define approximately, using Equation (5), as 



min[M 90 ]L n , k f 1 " e~ r( W ; r(t) < Rp 
e( ) " rn n c^ nJ X \ e -^ ; r(t) > Rp (6) 

The validity of this approximation depends on the fast rise and exponential decay of the function with time. 
Here tgoi n ,k is the total energy carried by the neutron outflow in case of a GRB central engine active for a 
time igo. Note that in case of complete decoupling of the neutron component from the rest of the outflow 
for F n j <C T p /, the observed GRB may result from the leading proton outflow with tgo duration and the 
lagging neutron component may have a different time scale. However, we ignore this possibility for the sake 
of simplicity. 
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The total energy released in a beta decay is the difference between the rest mass energies of the initial 
and final particles: Q = m n c 2 — (m p + m e + m u )c 2 = 0.782 — m v c 2 MeV [see, e.g., Krane (1987)]. This 
energy is shared by the kinetic energies of the electron (T e ), neutrino (T u ) and proton (T p ) in the decaying 
neutron's rest frame. With a negligible neutrino mass and tiny proton recoil energy: T p ~ 0.3 keV (because 
of its large mass), the maximum electron's kinetic energy is T e max ~ Q for T v — > 0. The shape of the 
kinetic energy spectrum of beta decay electrons is given by 

N e (T e ) = v / T 2 + 2T e m e c 2 (Q - T e ) 2 (T e + m e c 2 ) (7) 

which vanishes for T e = and Q. 

The kinetic energy of the electrons in Equation (6) is then distributed according to Equation (7) and 
we write a normalized kinetic energy distribution of beta decay electrons at an observed time t, in the jet 
outflow, as 

dK = K(t)K(Z) 

dT 'e f Q K(T>)dTi 

We verify this expression by numerically integrating the right hand side over T' e in the range 0-Q and 
comparing with Equation (6) for different t. Figure 1 shows the distribution in Equation (8) plotted for 
different parameters for long and short GRB jets at a time tp = Tp/2Y n j (assuming z ~ 0) when the total 
number of freshly created beta decay electrons in Equation (6) is the maximum. Note that M' e increases with 
£o, however, it is not directly proportional since T n j is not directly proportional to £ D [see Equation (2)]. 
We have listed different values of tp and T n j in Table 2 for different choice of r] and £ in case of long and 
short burst curves plotted in Figure 1 . 



4. Radiation spectrum and lightcurve 

Beta decay electrons are created suddenly with a peak kinetic energy of 0.25 MeV in the neutron's 
rest frame. These electrons may be thought to be created initially at rest and then accelerated to a constant 
high speed v = c/3 e within a short time interval r ~ hc/E e , following the uncertainty principle. Here 
E e = T e + m e c 2 is the total energy of the electron. Rapidly accelerating electrons then emit radiation (often 
called inner bremsstrahlung in the literature; see Petrosian & Ramaty (1972) for an early astrophysical 
application) similar to the bremsstrahlung. The radiation spectrum (total energy radiated per unit frequency 
interval) is also flat as bremsstrahlung and is given by (Jackson 1999) 



dl_ 

doj 



ire 



Pe 



1 + & 



(9) 



where the electron's speed (with momentum p e ) is 



iie ~ E e 



y/T e 2 + 2T e m e c 2 
T e + m e c 2 



0.74 



(10) 
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Fig. 1. — Kinetic energy distribution [Equation (8)] of beta decay electrons normalized to the total number 
of electrons [Equation (6)] in the GRB jet comoving frame at an observed time ta = Tp/2T n f for which 
r(t) = Rp or when the initial number of neutrons have decayed by a factor 1/e. The upper (lower) set of 
six curves correspond to long (short) GRBs. The solid, dashed and dotted-dash pair of curves correspond to 
7] = 100, 316 and 500 with £ G = 10 (£ G = 1) for the upper (lower) curve in each pair. The other parameters 
for both the long and short GRBs are listed in Table 1. The calculated values of T n j and tp for all curves 
are listed in Table 2. 
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The last number is for electrons of peak kinetic energy. Note that the spectrum defined in Equation (9) is 
proportional to the more familiar notation F v . 

The maximum energy of the emitted photons from beta decay electrons is huj max < E e following the 
condition lot < 1, although a full quantum mechanical treatment is necessary to evaluate the precise value. 
Semi-classically an electron loses a small fraction of its energy given by integrating Equation (9) over u as 



-E-y 

El. 



Q 



- 2 



0.001 



(11) 



Again, the last number is for an electron of peak kinetic energy 0.25 MeV. However, the collective radiation 
energy from all beta decay electrons in a GRB jet is large. For an observed GRB, the beta decay radiation is 
beamed along the line of sight of an observer because of the relativistic bulk motion of the jet even though 
beta decay electrons do not have any particular angular orientation and the radiation is emitted along the 
electron's directon in the jet comoving frame. 

We write the total energy radiated per unit frequency by all beta decay electrons at a time t in Equation 
(8) as would be measured by an observer as 



dX . s 



l + z)J dTi [t) dio' dle 



7rc(l + z) 



® dN' e 



dT: 



-In 



l + Pe 



- 2 



dT' 



(12) 



which is flat in frequency. The corresponding total radiation energy emitted over all frequencies at a time t, 
by integrating Equation (12) over < fou> < {T' e + m e c 2 )T n j/(\ + z), is 



q 2 T 2 



tot 



(t) ~ 



irhc(l + z) 



rQ 

/ {T> 
Jo 



l-P'e 



dT' e 



(13) 



We have plotted the total energy radiated by beta decay electrons in Figure 2 as function of observed 
time by numerically evaluating the integral in Equation (13) for different GRB model parameters used in 
Figure 1 and listed in Table 2. The total number of electrons at the peak kinetic energy plotted in Figure 
1 times the total energy radiated by each of these electrons using Equation (11) roughly corresponds to the 
peak total energy plotted in Figure 2 after multiplying with the corresponding T^ j factor. The maximum 
observed photon energy in all cases is e 7jmax ~ 1.3r n j/(l + z) MeV. 



5. Detection prospects and Discussion 

The highest energy photons carry most of the radiation energy plotted in Figure 2, since the beta decay 
radiation spectrum from Equation (12) is flat. The corresponding photon flux (energy received per unit area 
per unit time) on Earth, in a given energy band e 7i /-e 7jM ; is 



^( e j,l e 7,w) — 
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Table 1. Typical Gamma-Ray Burst Parameters 



Parameter 


Long GRB 


Short GRB 


Total isotropic-equivalent 


W 52 L 52 


10 50 L 50 


energy outflow (L) ergs/s 






Burst duration (tgo) s 


10 


1 


Redshift (z) 


1 


0.1 


Luminosity distance 21 (<i/J cm 


10 28 ^, 28 


10 27 (i Li27 


Baryon loading parameter (77) 


316r?2.5 


316??2.5 


Initial outflow radius (R a ) cm 




10 6 i? o ,6 



a Using typical cosmological parameters 




0.01 0.1 1 10 100 

Ms) 

Fig. 2. — Total radiation energy emitted by beta decay electrons over all frequencies as function of time 
[Equation (13)]. The model parameters for the GRB jets for different curves are the same as in Figure 1 and 
listed in Table 2. The upper (lower) set of curves correspond to the long (short) bursts. Note that the signal 
is stronger and peaks later for T n j = 77 cases as evidenced by the solid pair of curves in case of both long 
and short bursts. The dependence on £ is rather weak. 
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Upcoming high energy 7-ray detectors such as the Large Area Telescope (LAT) on board the Gamma- 
ray Large Area Telescope (GLAST) are sensitive to photons in the energy range 20 MeV - 300 GeV. The 
threshold fluence for LAT is ~ 4 x 10~ 8 ergs/cm 2 for a short integration time. This roughly corresponds to 
a few photons with energy > 0.1 GeV collected within a few seconds in its average 5000 cm 2 effective area 
over this energy range. We have listed in Table 2 the expected beta decay radiation flux at photon energies 
e 7 > 0.1 GeV, for different values of the relevant parameters. We used a luminosity distance cLl = 10 27 d,L,27 
cm corresponding to a redshift z ~ 0.1 for both the long and short bursts. This redshift is typical for short 
bursts, and is at the lower end for long bursts. (Observations indicate an occurrence rate of ~ 0.25 year -1 
for long bursts and 2-3 year -1 for short bursts in this redshift range). 

However, the beta decay signal can be diluted by other types of emission in the same GRB, since the 
signal can coincide temporally with the usual prompt 7-ray emission phase. This is somewhat mitigated 
by the fact that the two signals have different spectra, and the beta decay signal has a smoother rise and 
quick decay, compared to the sometimes erratic and longer signal of the long burst prompt emission. The 
GRB afterglow emission occurs on an even longer time scale than the beta decay time scale tp, and should 
not interfere with the latter. Hence we discuss the GRB prompt 7-ray emission as a source of background 
radiation below. 

The observed isotropic-equivalent bolometric 7-ray luminsosity, L 7 , of a long (short) GRB is ~ 
10 51 L 7 5i ergs/s (~ 10 49 L 7 4g ergs/s) which is mostly concentrated at the peak energy e 7iP k = 100e 7i _4 
keV. In the. fireball shock model, this energy corresponds to a fraction e e ~ L 7 /L ~ 0.1e e ,-i of the total en- 
ergy converted by shock accelerated electrons. The luminosity at higher e 7 decreases by a power-law, with 
index 2 < a < 3, following the synchrotron and/or inverse Compton (IC) radiation spectrum by a power- 
law distributed electrons. For a > 2, most of the energy in a given band e 7j j-e 7jtt is carried by the low 
energy photons (unlike in the beta decay signal). The corresponding prompt (background) synchrotron/IC 
photon flux at Earth is approximately 



for a long GRB at z ~ 0.1. Thus, the detection of the beta decay radiation signature would be background 
limited in the case of long bursts with £90 > 2 s, if the prompt spectrum extends to ~ 0.1 GeV with the above 
nominal spectrum. However, one can expect a sizable fraction of long GRBs where the prompt spectrum 
does not contain significant energy at 0.1 GeV, e.g. due to a spectral index which is steeper than the nominal 
value a = 3, or due to a fall-off in the spectrum [e.g., Preece, et al. (2000); also, as shown by Ryde (2004) 
that a fraction of long GRBs are dominated by a thermal-like component, with little or no emission above 
a few MeV]. While data in the ~ 0.1 GeV range is currently still limited, one may conservatively estimate 
at 30% level the fraction of all long GRBs which do not have substantial emission at 0.1 GeV, and hence 
which contribute no background for the beta decay signal. 

On the other hand, for short bursts, the beta decay signature peaks on timescales longer (see Figure 2) 




(15) 
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Table 2. Beta Decay Radiation Parameters 



Parameter T n j tp (s) JF a 



long GRB 



V = 


500 


£o = 


l 


376 


1.2 


0.33 






6> = 


10 


213 


2.1 


0.61 


V = 


316 


£o = 


1 


316 


1.4 


0.53 






& = 


10 


248 


1.8 


0.97 


V = 


100 


£0 = 


1 


100 


4.4 


1.7 






£0 = 


10 


100 


4.4 


3.1 



short GRB 



V = 


500 


£o = 


l 


174 


2.5 


2.6 






& = 


10 


99 


4.5 


1.4 


V = 


316 


£o = 


l 


203 


2.2 


2.4 








10 


115 


3.8 


4.9 


V = 


100 


£o = 


1 


100 


4.4 


3.8 






& = 


10 


100 


4.4 


6.9 



a in units of ( x 10 9 ) ergs cm 2 s 1 for long 
GRBs and in units of (xl0~ 12 ) ergs cm -2 s" 1 
for short GRBs 
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than the prompt emission (tgo < 2 s), and the beta decay signal is expected to be essentially free from a 
prompt emission background at 0. 1 GeV. 

One other source of background to consider is the diffuse extragalactic 7-ray background measured by 
EGRET (Sreekumar, et al. 1997). The corresponding photon number flux may be fitted in the 30 MeV-120 
GeV energy range as dN/de^ « 7.32 x 10~ 9 (e 7 /451 MeV) -21 cm" 2 s" 1 sr" 1 MeV" 1 . This background 
energy flux within ~ 1° angular resolution of LAT is F^ g ss 3 x 10~ 14 ergs cm -2 s" 1 for e 7 > 0.1 GeV 
which is negligible compared to the beta decay flux listed in Table 2 or GRB prompt flux in Equation (15). 

The predicted fluxes 0.1 GeV) in Table 2 were calculated for bursts of nominal (MeV) luminosi- 
ties Lj = 10 51 erg/s (long bursts) and L 7 = 10 49 erg/s (short bursts), at a redshift z = 0.1. Multiplying 
these nominal fluxes by a typical average beta decay peak duration tp ~ 5 s, and comparing to the LAT 
threshold sensitivity fluence 4 x 10~ 8 ergs/cm 2 at e 7 ~ 0.1 GeV, one sees that for short bursts, even though 
it is background free, the signal is undetectable by the LAT. For long bursts, however, it may be possible to 
detect the beta decay radiation signature from nearby long bursts whose luminosity is ten times the average 
value. Bursts satisfying these criteria may represent ~ 0.1-1% of a nominal total of 100 bursts per year 
detected by GLAST. This suggests that GLAST may be marginally able to detect such signals, and future 
generations of very energetic space-based 7-ray large area telescopes may be able to quantitatively explore 
this problem. 

Non-thermal multi GeV energy neutrino and 7-ray emission from the n-p decoupling phase can indicate 
the presence of a neutron component in the GRB jets, for values of the dimensionless entropy 77 > rj np , which 
is typically an unknown parameter. On the other hand, the electromagnetic radiation signature discussed here 
is valid for all values of 77. The evidence for a neutron component from the beta decay radiation signature 
would definitely imply the presence of protons in the GRB jet, since neutrons are coupled to them at least 
initially. Even though the converse is not true, the beta decay radiation signature may be an alternate way 
to explore the baryon loading in the GRB jets in the cases when high energy neutrinos from the internal 
or external shocks are absent (e.g., if the protons are not co-accelerated with electrons), or if too few are 
detected. The beta decay photon signature of neutron decay may thus be a valuable tool for investigating 
the particle acceleration process and possibly constrain the progenitors of Gamma-Ray Burst sources. 

Research supported in part through NSF AST 0307376 and NASA NAG5-13286. 
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